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Abstract. Due to the intrinsic properties of neutrinos, the gravitational lens effect 
for neutrino should be more colorful and meaningful than the normal lens effect of 
photon. Other than the oscillation experiments operated at terrestrial laboratory, 
■ in principle, we can propose a completely new astrophysical method to determine 

\ not only the nature of gravity and spacetime of lens objects but also the mixing 

£NJ ' parameters of neutrinos by analyzing neutrino trajectories near the central objects. 

However, compared with the contemporaneous telescopes through the observation of 
' the electromagnetic radiation, the angular, energy and time resolution of the neutrino 

telescopes are still comparatively poor, we just concentrate on the two classical tests 
5h , of general relativity, i.e. the angular deflection and time delay of neutrino by a lens 

object as a preparative work in this paper. In addition, some simple properties of 
neutrino lensing are investigated. 
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Neutrinos are one kind of the fundamental particles that make up the universe. Up 
to date, many interesting neutrino astrophysical sources have been known, such as the 
ultra-low energy Big Bang relic background neutrinos, the ~ 10 keV stellar neutrinos, 
the ~ 10 MeV supernova neutrinos, the ~ 100 TeV neutrinos from young supernova 
remnants in our galaxy, the ultra-high energy neutrinos from Active Galactic Nucles 
(AGNs), Gamma Ray Bursts (GRBs), Microquasars or some other celestial sources. 
Since neutrinos are electrically neutral leptons, they only participate in the weak 
interaction and are almost unaffected through transmission, and they open a new 
window into space to the high-energy processes in the universe to astronomy. As we 
known, several under-ice and under-water neutrino telescopes are being or scheduled 
to be constructed now, such as Icecube located at the south Pole worked in energy 
range from GeV to PeV, expected to achieve an angular resolution below 1°, a real-time 
resolution at 2 ns and an energy resolution below 0.3 in log(E). 

In addition, neutrinos are also hot topics and research areas for particle physicists. 
One of the most remarkable discoveries in particle physics connect with neutrinos in 
the last few decades is the finding that neutrinos are massive and mixed particles. 
Although the absolute scale of neutrino masses is still unknown, the evidence is strong. 
The precision measurement of neutrino masses and oscillation parameters now is one of 
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the main goals of the neutrino experiments. It is possible for researchers to measure the 

absolute mass of neutrino with an accuracy of 0.05 eV soon by the double beta decay 

experiment in many laboratories or by the lensing of the Cosmic Microwave Background 

radiation (CMB) from cosmological probes. 

As we known, a light ray can be deflected by gravity, it is same for neutrino. 

Since Soldner firstly published the derivation of the deflection of light by a massive 

object based on the framework of Newtonian gravity almost two centuries ago, ED the 

gravitational lensing has been full investigated numerically or analytically. Now, lensing 

of light ray has been used widely in the research of astrophysics and already became an 

excellently educational stuff for students of astronomy. ^ Another kind of gravitational 

lens researchers less mentioned is neutrino lensing, which has also been put forward and 

widely studied for a quite long time. Unlike photon, neutrino can penetrate into most of 

celestial objects easily except black holes or some special compact objects since it only 

participates in the gravitational and weak interaction. Therefore, one obvious advantage 

of neutrino lensing in principle is that it may provide a distinct method to probe the mass 

distribution of lens objects. Based on the standard solar model, Gerver investigated 

the possibility of neutrino focusing by the Sun's core. El Escribano et alM discussed 

the gravitational lensing of neutrinos by some extended sources as stars, galaxies and 

galactic halos. Moreover, Mena et alM studied the lensing effect of supernova neutrinos 

by the center black hole of our Galaxy. Eiroa and Romero^ analyzed the lensing effect 

of the cosmological sources of neutrinos by some supermassive black holes. However, 

all of works cited here used the lightray approximation, treated neutrino as a massless 

particle. Strictly speaking, neutrino does not move along null geodesies since it has a 
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non-zero rest mass. Based on this idea, Barrow and Subramanian LLI offered a possible 
explanation for the time delay of neutrinos from supernova SN1987a. 

In this Letter, we will show that this small but non-zero rest masses of neutrino will 
bring us some tinily small but non-trivial difference of the angular deflection and time 
delay of neutrino in gravitational field. In principle, neutrino lensing itself may provide 
a unique ground for probing the nature of neutrinos. Thus, we investigate gravitational 
lensing of neutrinos by the point-mass lens and re-investigate two famous effects of 
General Relativity, i.e. the deflection and the Shapiro time delay, we extend them to 
more general case. 

We start from the Kerr metric in Boyer-Lindquist (BL) coordinates 

ds 2 = p 2 ( — + d6 2 ) + (r 2 + a 2 ) sin 2 dd<f) 2 - c 2 dt 2 + -Ar (a sin 2 9d<j) - cdt) 2 , ( 1) 
A p 2 

where r s = 2GM/c 2 is the Schwarzschild radius, a = J/Mc is the specific angular 
momentum with the dimension of length, p 2 = r 2 + a 2 cos 2 9 and A = r 2 — r s r + a 2 . 
Since this metric describe a stationary and axisymmetric spacetime, it is easy for us to 
obtain those conserved quantities for test particles in Kerr geometry, 

E/mc = (1 - ^r)ci+ T -^ra sin 2 8<p , (2) 
p 2 p l 
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The dots denote differentiation with respect to an affine parameter. The constants E/m 
and L z /m can be viewed as the energy and the axial component of angular momentum 
per unit mass at infinity. In addition, other two constants of motion are the particle's 
rest mass m, 

mc = \p\ = {-g^VpPvf 12 , (4) 

and Carter's constant Q, 

Q = p* + cos 2 6 [a 2 (m 2 c 2 -E 2 /c 2 ) + (L/ sin 9) 2 ] . (5) 

For simplicity, we only investigate the motion for the open (or infinite) orbits in the 
equatorial plane, i.e. E > mc 2 , 9 = n/2 and pe = 0. Considering an incoming particles 
with Lorentz factor 7 = _E/mc 2 and impact parameter b L z j (V7 2 — Imc) at distance, 
we can solve for the velocity components in terms of 7 and b 
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It is obvious that the above equations are not only suitable for null geodesies for massless 
particles like the photon (7 — > 00) but also suitable for time-like geodesies for massive 
particles with finite 7. The orbit of a photon is more simple since it is not decided by its 
energy but completely decided by its impact parameter. We will focus on the general 
behavior of test particles (including photons) hereafter. 

Setting dr/d(j) = 0, we can derive a relationship between the turning point (the 
point of closest approach to central object) r and the impact parameter b, 
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where + and — denote prograde and retrograde particles respectively. We also can 
express r as a power series in u — r s /\b\, 
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Figure 1. Deflection angle for photon or some other relativistic particles with different 
impact parameters. The solid and the dashed curves denote the numerical solutions 
and the series solutions, respectively. The circle and the triangle curves represent the 
deflection of the retrograde and the prograde particles (with s = 1 as a demonstration) 
in the equatorial plane of the extreme Kerr Black Hole. 



where s = ^7 2 — 1? 5 = a/r s is the specific angular momentum in units of the 
Schwarzschild radius, with a > for prograde particles and a < for retrograde particles 
in a unified expression. 

Simply and repeatedly used the mathematical relation d u fo^g(x,u)dx = 

lo^ d u [g(x, u)] dx + f'(u) g(f(u),u), the deflection angle (in radians) can be expressed 
as a Taylor series about u as 
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If you do not care about the interminable derivation and integration, definitely you can 
expand the above integration to any order you want. As a demonstration, geodesies are 
expanded as a Taylor series up to and including fourth-order terms in r s /\b\ here. When 
a = and s — > oo, we recover the famous first order light deflection in Schwarzschild 
space-time, 2r s /b, and add several high order modification. Fig.l shows the logarithm 
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plot of the deflection angle for photon and some other relativistic particles with different 
impact parameters. For the convenience of comparison, we have carried out the 
deflection integral numerically, i.e. the solid curves in Fig.l. Clearly, the Taylor series 
expansion offered here are quite well formula to specify the bending angle especially 
for small angle, weak field (with \b\ much larger than r s ) and ultra-relativistic particles 
(with s ^> 1). The intermediate relativistic particles with s — 1 are taken as an example 
to demonstrate the rotational contribution. It is natural and obvious that the spin is 
a high-order modification, which only appears in second-order or even higher order 
terms of the expansion, which means that the spin parameter of black hole only can be 
determined by the careful observation of the large angle scattering of the intermediate 
or ultra relativistic particles. Due to the frame-dragging effect of the rotating black 
holes, the particle traveling in the direction of rotation of the object will move around 
the central object faster than particles moving against the rotation, the position of the 
turning point for prograde particle will be even more far away from the center than 
that of retrograde particles. Thereupon the deflection angle is enhanced for retrograde 
particles (the circle curves) and reduced for prograde particles (the triangle curves). 
The vertical dashed line in the infixed plot, r s /b = 2/y27, corresponding to the critical 
impact parameters b c for Schwarzschild black hole. For b < b c , the photon will spiral 
in and be captured by the central black hole. For more general case in Schwarzschild 
spacetime, the critical value is 

b _ r VW - 2) 2 (3 7 2 - 4) + 7v /9^8 (9 7 4 - 20 7 2 + 12) 
v / 2( 7 2 -l)( v / 9 7 2 - 8 - 7 ) 

However there is no analytic expression for the critical impact parameter in Kerr black 
hole. By considering the capture of the background photons, hot, warm and cold dark 
matter in the early Universe, we plan to investigate the growth of primordial black holes 
in our future work, which are possible seeds for supermassive black holes at the center 
of the galaxies. 

We have discussed the angular deflection of photon or relativistic particle by a 
rotating spherical body. From the first order term of Eq. (11), we can see that the 
difference of the deflection between photon and ultra-relativistic particle, 2r s /6( 7 2 — 1), 
is tinily small, which means that the neutrino lens and the photon lens are nearly 
the same in weak deflection limit, that is why people can safely and simply use null 
geodesies as photons do to study neutrino lensing.^ ^ If we want to directly determine 
the mass of neutrinos by bending, we must seek for the observation of the large angle 
scattering. Unlike photons, neutrinos are electrical neutrality, only participate in the 
weak interaction, can penetrate the lens objects easily, hence it should be possible for 
researchers to detect those large angle scattering neutrinos in principle. However, it is 
difficult for us to find the optical counterpart of the same source to compare with due 
to the strong dust extinction in the neighborhood of the lens objects. Furthermore, 
the requirement for angular resolution may go far beyond the capability of neutrino 
telescopes at present even in the near future, the operation process to detect the bending 
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Figure 2. Sketch of a typical lensing system. 



angle precisely will be very difficult and nearly hopeless for us. Thus, it seems that our 
work here is more meaningful from the theoretical than from the observational point 
of view. In fact, we really plan to study the neutrino lensing by a Kerr black hole in 
weak deflection limit with higher order modification by following Sereno's wor m 
a forthcoming paper. In the following paragraphs, we only devote to the study of the 
properties of neutrino lensing as the simplest, i.e. the point-mass lens model is employed. 
In addition, we only consider the contribution of the first order term for simplicity. 

If there are no other deflectors except a lensing object of mass M close to the 
line-of-sight to a source, the sketch map of a typical lensing system can be drawn 
easily (See Fig. 2). The distances between source and observer, source and lens, and 
observer and lens are given by Dos, Dls, Dol, respectively. The angle (3 and 9 
denote the positions of source and images with respect to the lens direction. If the 
physical size of the lensing object is much smaller than both Dls an d Dql, and the 
all angles (f3, 9, 5(f)) are very small, from Fig. 2 we can reach the geometric relation 
[3 = 9 — 5(f) D L s/ Dps- Notice b = 9D OL and the deflection angle expression to first 
order, we obtain the famous lens equation (3 = 9 — 9^/9, where the Einstein angle 
9e now satisfies 9\ = (2 + l/s 2 ) r s D L s/(DosDoL), which defines the angular scale 
for a lens situation. By solving the lens equation, one can find the positions of the 
primary (upper sign) and the secondary (lower sign) images for an isolated point source, 



9i,2 = \ \J3 ± y (3 2 + 4:9^ J , which correspond to the shortest geodesic and the longest 

geodesic with b\ = 9\Dol and 62 = \9-2\Dol, respectively. 

As we known, times can be measured to much greater accuracy than angles, we will 
further investigate another prominent and equivalent effect of General Relativity, i.e. 
the Shapiro time delay or gravitational time delay here. This effect was first pointed by 
Shapiro in 1964 about the time delay of signals passing near a massive object.^ Let us 
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calculate the time required for a particle to go form a source point with r = r$, 4> = 4>s, 
to an observer at r = ro, 4> — 4>o m the equatorial plane. The equation governing the 
time history of orbits is easy to be given by Eq. (6) and (8). The time required for 
particle to go from the turning point r to r is 



t(r, r ) 




r 2 r 2 + a 2 (2r 2 — Tq) — 2ar^J (r 2 — r%) (a 2 + Tq) 

We only keep the first order in r s /r in the integrand function. The result show that 
the traveled the time has two different components. The leading term is related 
to the geometrical of the spacetime in the absence of the lens, now is the time 
for the relativistic particle traveled in straight lines at unit velocity, c^'-f 2 — I/7, in 
the Minkowski spacetime. The other terms are the first-order approximation for the 
gravitational time delay. This is the well-known 'Shapiro effect', which has been amply 
tested by radar echo delay experiments in our Solar System. Of course, the total time 
required for a particle to go form a source point to an observer can be written as 



tos = t(D O L,r ) + t(yjDl s + Dq S (3 2 , r ). 

Since particles that form distinct images are almost emitted at the same time from 
the same source but travel by different paths, they may reach an observer at different 
time. Combining Eq. (13) with the above image equation, we can infer the physical time 
delay function between the gravitationally lensed images, 
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where f3 = (3/9e is the reduced misalignment angle. Furthermore, we can give the time 
delay of a relativistic particle relative to a photon emitted by the same source at the 
same time, At p h, which can be measured very easy. As space is limited, we shall not go 
into details of the formula. 

Besides the multiplicity and distortion of images and the time delay between 
different images, another prominent effect of lensing is the magnifications of the images, 
which are given by the ratio of the solid angles of the image and the source, 
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The measurable total magnification fi = fi + + //_ is always larger than one. Since the 
Einstein angle 9e now is a slightly larger than that of photon, the image magnification 
of neutrino lensing is slightly larger than that of normal lensing. As usual, the maximal 
magnification is infinite as (3 — > 0, but this is not a real situation since the physical 
objects should have a finite size R*, the effective limit now is /? — >• = R^/ Dos- 

Based on 17- year radial velocity and 12-year astrometric measurements of the short 
period star SO-2, last year Ghez et a/.^3 reported the newest and the most precise results 
about the central supermassive black hole in our galaxy. The balk hole is 8.0 ± 0.6 kpc 
from us with mass (4.1 ±0.6) x 1O 6 M . Supposed a type of two supernova explode at 8 
kpc away from our galactical center with anggular position (3 = 0.1 ~ 0.5 6e- The order 
of magnitude of the time delay between the primary and secondary images (no matter 
for photons or neutrinos) is several dozens seconds. Recently experimental results (from 
KamLAND, MINOS, Super-Kamiokande, etc.) and CMB observation show neutrinos 
have sub-eV mass (at least one mass eigenstate state neutrino with a mass of at least 
0.04 eV, the best estimate of the mass square difference are Am^ = 7.684 x 10~ 5 eV 2 , 
Am 2 32 = 0.0027eV 2 ).M-[l2],[l3] Taken E v ~ lOMev, the typical value of supernova 
neutrinos, m u c 2 ~ 1 eV, the flight time delay of neutrino relative to photon is about 
~ 10~ 5 s x^kpc \f^ r s ^ n2 @- Thus, we can really use the forthcoming neutrino facilities 
to do some more detailed analysis of time-varying. Combined with the observation of 
the normal lensing of light ray, neutrino lensing probably can provide us some valuable 
hints about the spacetime characteristic of the lens object, even the intrinsic quality 
of neutrino itself. However we only investigate the point-mass lens model, while the 
matter distribution during the propagation may has a strong impact on the results, so 
a more detailed and sophisticated treatments are need in the future. 
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